Pressure-induced Phonon Softenings and the Structural and Magnetic Transitions in 
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To investigate the pressure- induced structural transitions of chromium dioxide (Cr02), phonon 
dispersions and total energy band structures are calculated as a function of pressure. The first 
structural transition has been confirmed at P~ 10 GPa from the ground state tetragonal CrC>2 
(t-CrCh) of rutile type to orthorhombic Cr02 (o-Cr02) of CaCl2 type. The half-metallic property 
is found to be preserved in o-Cr02. The softening of Raman-active Bi s phonon mode, which is 
responsible for this structural transition, is demonstrated. The second structural transition is found 
to occur for P> 61.1 GPa from ferromagnetic (FM) o-Cr02 to nonmagnetic (NM) monoclinic Cr02 
(m-CrO^) of M0O2 type, which is related to the softening mode at q = R(i, 0,|). The third 



structural transition has been newly identified at P= 
type that is a FM insulator. 

PACS numbers: 61.50.Ks, 63.20.dk, 71.15.Rf, 74.25.Kc 



GPa from m-Cr02 to cubic Cr02 of CaF2 



I. INTRODUCTION 

Chromium dioxide (CrC^), which crystallizes in the 
tetragonal structure of rutile-type, is a well-known mate- 
rial because of its half- metallic nature with T c — 390 K. 1 
The origin of the ferromagnetic (FM) and half-metallic 
property of Cr02 was explained in terms of the double- 
exchange model3— Due to the crystal field of distorted 
(flattened) Cr-C>6 octahedra, Cr t2 g states are split into 
lower d xy and higher d xz /d yz states. Out of two d elec- 
trons of Cr 4+ , one occupies the lower d xy that is local- 
ized, while the other occupies the higher d xz /d yz that 
are delocalized near the Fermi level (Ep) due to the hy- 
bridization with O p states. Then the double-exchange 
interaction arises from the Hund coupling between local- 
ized d xy and delocalized half- filled d xz /d yz states, so as 
to produce the FM and half-metallic properties. 

In contrast to numerous reports on electronic and mag- 
netic properties of Cr02, there have been relatively small 
number of studies on structural and lattice dynamical 
properties of Cr02- Especially, there is no experimen- 
tal or theoretical report on the phonon dispersion curve 
for Cr02, except for a few Raman studies4r£ Under the 
pressure, Cr02 is known to undergo the structural tran- 
sition from the ground state tetragonal Cr02 (t-Cr02) 
to the orthorhombic &O2 (o-CrC^) of CaCki-type at P 
= 12-14 Gpa^ii The question followed is whether there 
will an additional structural transitions from o-Cr02 at 
higher pressure. In fact, this question is not just for 
CrC>2 but also relevant to the structural stability issue 
of transition-metal (TM) dioxides (TM0 2 ). Note that 
TMO2S show diverse structures depending on the TM 
elements^— Furthermore, magnetic properties of &O2 
under the pressure are intriguing, such as (i) whether 
the half-metallic nature is preserved, and (ii) when &O2 
becomes non-magnetic. 

In this work, to investigate the pressure-induced struc- 
tural transitions of Cr02, we have studied phonon disper- 
sions and total energies of relevant &O2 structures as a 



function of pressure. Based on the calculated phonon 
dispersions and the total energies, we have found three 
possible structural transitions with increasing pressure. 
The first transition is consistent with the known transi- 
tion from t-CrC-2 to o-CrC>2. At this transition, FM and 
half-metallic properties are preserved, in agreement with 
previous reports of literatures^ The second transition 
is from o-Cr02 to monoclinic Cr02 (m-CrC^) of M0O2 
type, which is nonmagnetic (NM). The third transition is 
identified from m-Cr02 to cubic CrC>2 (c-CrC^) of CaF2- 
type. Interestingly, c-Cr02 is a FM insulator even at the 
high pressure of P > 88.8 GPa. Note that the second 
and third structural transitions are our new findings for 
Cr02 under the high pressure. 



II. COMPUTATIONAL DETAILS 

Band structures and phonon dispersions of Cr02 
were obtained by employing the pseudo-potential band 
method and the linear response method, respectively, im- 
plemented in the Quantum ESPRESSO code.i^ The 
generalized gradient approximation (GGA) is utilized for 
the exchange-correlation potential. Self-consistent elec- 
tron and phonon band calculations were carried out after 
the full-relaxation of internal atomic positions and lattice 
parameters. 

We have considered various structures of Cr02- At the 
ambient pressure, the stable phase is t-Cr02 of rutile- 
type (P^/mnm), in which Cr atoms are positioned 
at (0,0,0) and (|, |, |), while O atoms at ±(u, u, 
0) and ±(i+u, |— u, |). Initial lattice constants and 
atomic positions adopted before the full-relaxation are 
a=b=4.42lA, c=2.916A, and u=0.3043^ For the high 
pressure phase of o-Cr02 of CaCi2-type (Pnnm), we have 
adopted a=4.3874A, b=4.2818A, c=2.8779l, u x =0.299, 
and u y =0.272j^ For candidate structural phases at the 
higher pressure, we considered m-Cr02 of M0O2 type 
(P2i/c)±&±± and c-Cr0 2 of CaF 2 -type (Fmlm)^ In the 
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lattice displacements induce the structural transforma- 
tion from t-CrC>2 of rutile type to o-Cr02 of CaCl2 type. 
Figure [T]Jc) provides the phonon dispersion of o-CrC>2 at 
P=34.1 GPa. The phonon dispersion is regular, implying 
that o-CrC>2 is stable at this pressure. Therefore, Fig. Q] 
clearly demonstrates that there is a structural transition 
from t-CrC>2 to o-CrC>2 at the pressure of P< 34 GPa. 

In Fig. EJa), we plotted the calculated Raman-active 
phonons of t-Cr02 and o-Cr02 as a function of pressure. 
There are four Raman modes (Bi g , E g , Ai g , B2 g ) for 
t-CrC>2, and six Raman modes (A g , Bi g , B2 g , B3 g , A g , 
Bi g ) for o-Cr02i^ ' 16 ' 17 Our data are consistent with ex- 
perimental data up to Prj 40 Gpa; 4 With increasing the 
pressure, one can clearly see the softening of Bi g mode 
of t-Cr02, which indicates the structural instability of 
t-Cr02- One can also notice two transition points. The 
first one corresponds to the transition from t-Cr02 to 
o-Cr02 at P=9.8 GPa. At this transition, Cr02 keeps 
its FM and half-metallic properties 4iLIii. The second 
one corresponds to the transition at P=76.0 GPa. The 
stable structure for P> 76 has not been identified yet. 

The phonon anomalies at two transition points are also 
revealed in the variation of lattice constants of Cr02 un- 
der the pressure. In Fig. [3fb), the lattice constants are 
plotted as a function of pressure. One can see anomalous 
behaviors of the lattice constants at the two transition 
points that are coincident with those in Fig. [3{a). 
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FIG. 1: (Color online) The phonon dispersion curves of Cr02- 
(a) FM t-Cr02 of rutile type at the ambient pressure, (b) FM 
t-Cr02 at P=34.0 GPa. Notice the phonon softening at T, 
which corresponds to the Bi 9 mode. The negative frequency 
here represents the imaginary part of the phonon frequency, 
(c) FM o-CrQ 2 of CaCk type at P=34.1 GPa. 



latter, a Cr atom is positioned at (0, 0, 0), and O atoms 
at (0.25, 0.25, 0.25), (0.25, 0.25, 0.75). 



III. RESULTS 



Figure Q] shows the phonon dispersions of t-Cr02 and 
o-Cr02 at the ambient and high pressures. As shown 
in Fig. [lja), t-Cr02 at the ambient pressure has regular 
phonon dispersions, reflecting the stable phase of t-Cr02 
at the ambient pressure. In contrast, t-Cr02 at P=34.0 
GPa in Fig.[IJb) has a softening phonon mode at q=r, in- 
dicating the structural instability of t-Cr02 at this pres- 
sure. The softened mode corresponds to B lg mode that 
is Raman- active. As shown in Fig. [5J Bi g mode gener- 
ates the rotating motions of oxygen ions. The resulting 



To investigate the second structural transition in more 
detail, we have examined the behavior of magnetic mo- 
ment. Srivastava et al^ once reported that there would 
be a magnetic transition in t-Cr02 from half-metallic to 
NM at P« 65 GPa. However, as discussed above, there 
occurs a structural transition from t-Cr02 to o-Cr02 at 
the low pressure of about P=10 GPa. Hence, in Fig.^Ja), 
we have examined the magnetic moment behavior for o- 
Cr02- It is seen that there is a FM to NM transition 
at P=76.0 GPa, which is close to the second structural 
transition point. 
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FIG. 2: (Color online) The normal mode of the Bi 9 soft 
phonon at F for P=34.0 GPa. The blue and light-blue circles 
represent Cr ions, while the orange and red circles represent 
oxygen ions. Blue and orange ions are located at z=|. Only 
the oxygen ions move in this mode. 
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FIG. 3: (Color online) (a) Calculated Raman-active phonon 
frequencies of t-Cr02 and o-Cr02 versus pressure, t- and o- 
stand for t-Cr02 and o-Cr02 respectively. The lines connect- 
ing data are guide for eyes. Two phase transitions are noticed 
at P=9.8 GPa and P=76.0 GPa, which are marked by vertical 
lines, (b) Calculated equilibrium lattice constants of t-Cr02 
and o-Cr02 versus pressure. 



The magnetic moment of o-CrC>2 suddenly drops at 
this transition point. The half- metallic property persists 
up to this pressure. This magnetic transition was also 
observed by Kuznetsov et alZ, who obtained the transi- 
tion pressure of P=53 GPa based on the pseudo-potential 
band method implemented in the VASP code. 

It is thus tempting to conjecture that the second tran- 
sition observed in Fig. EJa) corresponds to the magnetic 
transition in o-Cr02- However, the phonon dispersion 
curve in Fig. gjb) for NM o-Cr0 2 at P=80.5 GPa shows 
the strong phonon softenings, indicating that even NM o- 
Cr02 is unstable at the pressure of P> 76.0 GPa. There- 
fore, it is not possible that FM o-Cr02 changes into NM 
o-Cr02 with increasing the pressure. There might be an 
additional structural transition in this pressure range. 



Two candidate structures after the transition are m- 
Cr0 2 of Mo0 2 type and c-Cr0 2 of CaF 2 type. The 
monoclinic structure of Mo0 2 type is chosen from the 
expectation that, with increasing pressure, two d elec- 
trons of Cr become itinerant, and the local environment 
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FIG. 4: (Color online) (a) Calculated magnetic moment of 
o-Cr02 versus pressure. Magnetic transition from the FM 
to the NM phase occurs at P=76.0 GPa. (b) The phonon 
dispersion curve of NM o-Cr0 2 at P=80.5 GPa. 



becomes similar to that of Mo0 2 . In a similar system 
V0 2 , a softening of phonon frequency was observed at 
q=R (|, 0, i), which corresponds to the atomic move- 
ments from the tetragonal structure of rutile-type to the 
monoclinic structure.- 8 Indeed, the softening mode at 
q = R in Fig. SJb) is related to this structural transi- 
tion, because the orthorhombic structure of CaCl2 type is 
nothing but the distorted rutile-type structure. Figure [5] 
depicts the normal mode of the softened phonon mode 
at q = R. The displacements generate Cr-Cr dimeriza- 
tion along the c-axis, which is consistent with the main 
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FIG. 5: (Color online) The normal mode of the softened 
phonon at q = R for P=80.5 GPa. The blue and light-blue 
circles represent Cr ions, while the orange and red circles rep- 
resent oxygen ions. 
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FIG. 6: (Color online) Total energies of various Cr02 struc- 
tures versus volume. Data are fitted by the Birch-Murnaghan 
equation of state. 



to those in t-CrC>2 and o-Cr02- Note, however, that c- 
CrC>2 is an insulator not a half-metal. We have confirmed 
this result by employing the all-electron full potential lin- 
earized augmented plane wave method (FLAPW) band 
method 21 implemented in WIEN2k package too^ The 
present result is different from that by Srivastava et al.^ 
who obtained the stable NM metallic phase of c-Cr02 for 
P>90 GPa. The different result is likely to come from 
their use of a simple tight-binding LMTO band method. 

The transition from m-Cr02 to c-CrC>2 is thought to 
originate from the increasing packing ratio. There are 
six and eight oxygens around Cr in m-Cr02 and c-CrC>2, 
respectively. Haines et al~2- proposed several possible 
paths of structural transition from rutile to CaF2 type 
structure in TMC>2- Interestingly, one of the paths is the 
same as the present structural transition path, rutile- 
type {P^il ' mnm) — > CaCl2-type (Pnnm) —> Mo02-type 
(P2i/c) — > CaF 2 -type (Fm3m). But they did not take 
into account the magnetic state. 



distortions of the transition from the rutile-type to the 
monoclinic structure. It is thus reasonable to expect the 
second transition to be from o-Cr02 of CaCl2 type to 
m-Cr02 of M0O2 type. Concerning another candidate, 
c-Cr02 of CaF2 type, there were previous reports predict- 
ing the structural transitions from CaCl2 type to CaF2 
type structure for CrO^ and Ru02 - - Also, CaF 2 type 
is a typical structure of TM02- For example, ZrC>2 and 
Hf02, which have relatively large cations, crystallize in 
CaF2 type structure at high temperature^ 



To identify the additional structural transition at the 
higher pressure, we have compared the total energies of 
candidate structures in Fig. [6] From the total energy ver- 
sus volume curves in Fig. [51 one can identify three struc- 
tural phase transitions. The first one is from FM t-CrC>2 
to FM o-Cr02 at the estimated pressure of P=12.2 GPa, 
which is consistent with phonon calculation in Fig. [3] For 
P > 61.1 GPa, NM m-Cr02 becomes the most stable, 
which corresponds to the second transition from o-Cr02 
to m-Cr02, as discussed in the phonon study of Fig. S) 
For P > 88.8 GPa, c-Cr0 2 , which is a FM insulator, 
becomes the most stable. The more stable FM and in- 
sulating phase than the NM metallic phase of c-Cr02 at 
this high pressure is extraordinary. The magnetic mo- 
ment of c-Cr02 amounts to ~ 2/is/Cr, which is close 



IV. CONCLUSION 

We have studied the pressure effect on the structural 
properties of Cr02 by performing the phonon dispersion 
and total energy band structure calculations. Combin- 
ing two analysis methods, we have found that there are 
three structural transitions with increasing pressure up 
to 100 GPa. The first one is the transition from t-Cr02 
of rutile type to o-Cr02 of CaCl2 type at Pss 10 GPa 
(9.8 GPa from phonon dispersion analysis, while 12.2 
GPa from total energy study) . The FM and half- metallic 
properties of Cr02 are preserved at this transition. The 
second structural transition is from FM o-Cr02 to NM 
m-Cr02, which corresponds to the lattice displacement 
of the phonon softening at R in o-Cr02- The transition 
pressure is P=76.0 GPa from the phonon dispersion anal- 
ysis, whereas P=61.1 GPa from the total energy study. 
The third structural transition is from NM m-Cr02 to 
FM c-Cr02 at P= 88.8 GPa, which is accompanied by 
the metal to insulator transition. 



Acknowledgments 

This work was supported by the NRF (No. 2009- 
0079947), and the KISTI supercomputing center (No. 
KSC-2011-C2-36). 



1 K. Schwarz, J. Phys. F: Met. Phys. 16, L211 (1986). 

2 M. A. Korotin, V. I. Anisimov, D. I. Khomskii, and G. A. 
Sawatzky, Phys. Rev. Lett. 80, 4305 (1998). 

3 M. Katsnelson, V. Irkhin, L. Chioncel, A. Lichtenstein, 
and R. de Groot, Rev. Mod. Phys. 80, 315 (2008). 

4 B. R. Maddox, C. S. Yoo, D. Kasinathan, W. E. Pickett, 



and R. T. Scalettar, Phys. Rev. B 73, 144111 (2006). 

5 T. Yu, Z. X. Shen, W. X. Sun, J. Y. Lin, and J. Ding, J. 
Phys.: Condens. Matter 15, L213 (2003). 

6 M. N. Iliev, A. P. Litvinchuk, H.-G. Lee, C. W. Chu, A. 
Barry, and J. M. D. Coey, Phys. Rev. B 60, 33 (1999). 

7 A. Y. Kuznctsov, J. S. de Almeida, L. Dubrovinsky, R. 



5 



Ahuja, S. K. Kwon, I. Kantor, A. Kantor, and N. Guignot, 
J. Appl. Phys. 99, 053909 (2006). 

8 L. Mattheiss, Phys. Rev. B 13, 2433 (1976). 

9 R. Pynn, J. Axe, and R. Thomas, Phys. Rev. B 13, 2965 
(1976). 

V. Eyert, R. Horny, K. Hock, and S. Horn, J. Phys.: Con- 
dens. Matter 12, 4923 (2000). 

1 V. Eyert, Ann. Phys. (Leipzig) 11, 650 (2002). 

2 Quantum ESPRESSO (opEn-Source Package for Research 
in Electronic Structure, Simulation, and Optimization). P. 
Giannozzi et al., J. Phys.: Condens. Matter 21, 395502 
(2009). http://www.quantum-espresso.org 

3 We have used the kinetic energy cutoff for wave functions 
of ~ 400 eV. We have selected k-point samplings such that 
8x8x8 for t-Cr0 2 and o-Cr0 2 , 7x7x5 for m-Cr0 2 , and 
12x12x8 for c-Cr0 2 in the Monkhorst-Pack grid to set 
similar k-point density in the Brillouin zone of each struc- 
ture. 



W. Xue-Wei, N. Dong-Lin, and L. Xiao-Jun, Chinese Phys. 
Lett. 24, 3509 (2007). 

V. Srivastava, M. Rajagopalan, and S.P. Sanyal, The Eu- 
ropean Phys. J. B 61, 131 (2008). 

16 S. S. Rosenblum, W. H. Weber, and B. L. Chamberland, 
Phys. Rev. B 56, 529 (1997). 

17 W. H. Weber, G. W. Graham, and J. R. McBride, Phys. 
Rev. B 42, 10969 (1990). 

18 F. Gervais and W. Kress, Phys. Rev. B 31, 4809 (1985). 

19 J. Haines and J. M. Leger, Phys. Rev. B 48, 13344 (1993). 

20 A. Munoz-Paez, J. Chem. Edu. 71, 381 (1994). 

21 B. I. Min, H. J. F. Jansen, and A. J. Freeman, Phys. Rev. 
B 33, 6383 (1986). 

22 B. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, and 
J. Luitz, WIEN2K (Technische Universitat Wien, Austria, 
2001). 



